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Abstract. We present an investigation of galaxy-galaxy interactions and their eflects on the velocity fields of disc 
galaxies in combined N-body/hydrodynamic simulations, which include cooling, star formation with feedback, and 
galactic winds. Rotation curves (RCs) of the gas are extracted from these simulations in a way that follows the 
procedure applied to observations of distant, small, and faint galaxies as closely as possible. We show that galaxy- 
galaxy mergers and fly-bys disturb the velocity fields significantly and hence the RCs of the interacting galaxies, 
leading to asymmetries and distortions in the RCs. Typical features of disturbed kinematics are significantly rising 
or falling profiles in the direction of the companion galaxy and pronounced bumps in the RCs. In addition, tidal 
tails can leave strong imprints on the rotation curve. All these features are observable for intermediate redshift 
galaxies, on which we focus our investigations. We use a quantitative measure for the asymmetry of rotation curves 
to show that the appearance of these distortions strongly depends on the viewing angle. We also find in this way 
that the velocity fields settle back into relatively undisturbed equilibrium states after unequal mass mergers and 
fiy-bys. About 1 Gyr after the first encounter, the RCs show no severe distortions anymore. These results are 
consistent with previous theoretical and observational studies. As an illustration of our results, we compare our 
simulated velocity fields and direct images with rotation curves from VLT/FORS spectroscopy and ACS images 
of a cluster at z=0.53 and find remarkable similarities. 
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1. Introduction 

The evolution of galaxies within the dense environment of 
a galaxy cluster is strongly influenced by various interac- 
tion phenomena. These environmental effects are super- 
posed on the hierarchical growth predicted by cold dark 
matter (CDM) theories. Tidal interactions, for instance, 
have a higher probability in galaxy clusters and groups, 
due to the higher number density of galaxies. Going to 
higher redshifts, the merger frequency increases (Dressier 
et al. 1994), as galaxy clusters are assumed to assemble 
rather late in hierarchical CDM models. The relative ve- 
locities during encounters in galaxy clusters are higher 
than in the field and therefore reduce the probability for 
mergers. 

In a recent paper we studied the influence of tidal in- 
teractions and mergers on the star formation rates of the 
involved galaxies using numerical simulations (Kapferer et 
al. 2005). However, not only is the star formation affected 
by galaxy-galaxy mergers, but also the internal kinematics 
of the gas. Most of the irregularities in the velocity held 
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reflect disturbances in the gravitational potential of the 
galaxy. It is expected that kinematic disturbances fade 
within a few rotation cycles (< 1 Gyr) (e.g. Dale et al. 
2001). Therefore, they can be used to trace the recent in- 
teraction history of a galaxy and give possible clues to the 
type of the respective interaction. In addition, identify- 
ing disturbances in the RCs is important for TuUy-Fisher 
studies in order to accurately derive the maximum rota- 
tion velocity V^ax- 

Observations by Rubin, Waterman, & Kenney (1999) 
have shown, for example, that rotation curves from sev- 
eral Virgo cluster galaxies exhibit clear kinematic distur- 
bances. However, not all peculiarities are severe. Some 
rotation curves show only small asymmetries and others 
have no distortions at all. The authors did not flnd any 
strong correlations with Hubble type, local galaxy density, 
or HI deficiency. Many other studies have also relied on 
ID information about the internal kinematics and used 
RCs to trace interactions in galaxy groups and clusters 
(e.g. Swaters et al. 1999; Dale et al. 2001; Conselice & 
Gallagher 1999; Mendes de Olivera et al. 2003). During 
the past few years it has become feasible to observe the 
full 2D velocity field of local galaxies. For example Chemin 
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et al. (2006) investigated the Ha line of 30 spiral and ir- 
regular galaxies in the Virgo cluster using Fabry-Perot 
interferometry (FPI) . Their images show different signs of 
interactions in several objects. Likewise, the GHASP sur- 
vey observed local field galaxies (Garrido et al. l20U2|l with 
FPL 

Almost all studies of distant, faint, and small galaxies 
are still based on slit spectroscopy (Vogt 2001, Bohm et al. 
2004). Focussing on the environment of clusters, Ziegler et 
al. (2003) observed spirals in three intermediate redshift 
clusters (z=0.3-0.6), using multi-object spectroscopy with 
VLT/FORS. In the same manner, Bamford et al. (2005) 
investigated five distant clusters. Using Subaru/FOCAS, 
Nakamura et al. (2006) examined four distant clusters. 
Finally, Metevier et al. (2006) studied a cluster at z=0.4 
using Keck/LRIS. All these studies analysed the luminos- 
ity evolution of the galaxies via the Tully-Fisher relation. 
While Ziegler et al. (2003) and Nakamura et al. (2006) do 
not find significant differences between galaxies in clusters 
and in the field, Bamford et al. (2005) claim that cluster 
galaxies are on average brighter than their field counter- 
parts. Part of this discrepancy may be attributed to the 
way distortions and irregularities in the RCs are accounted 
for. 

Therefore, to systematically explore the influence of 
interactions on the velocity fields and hence the RCs, 
we present here a comprehensive analysis using numeri- 
cal simulations. We present the results by closely mim- 
icking the observations cited above. To this purpose we 
extract RCs that can be directly compared to the ob- 
served ones. Our N-body/SPH simulations have already 
been successfully applied to model RCs of isolated spiral 
galaxies (Kapferer et al., 2006). In that work we exam- 
ined, for example, how observational constraints affect the 
quality of RCs. Here we go one step further to investigate 
the kinematics of tidally interacting and merging systems. 
Depending on the interaction geometry, galaxy mergers 
and fly-bys will cause distortions of different strengths in 
the velocity field of galaxies. 

The main questions we will answer are what kind of 
distortions in the RCs can arise from tidal interactions? 
Are these peculiarities observable at intermediate and high 
redshifts? Are the features permanently visible or do they 
vanish within certain timescales? Some of the questions 
that we investigate in the present work were also ad- 
dressed by Barton et al. I|1HHH|I and Heyl et al. (1996), 
who found qualitatively similar results to ours. However, 
as both groups used purely N-body models, they were not 
able to extract RCs as close to the observational procedure 
as we do (see Sect. 2.1). This is especially crucial in the 
treatment of high redshift galaxies, which we focus on in 
the present study. A very interesting work connecting ob- 
servations and simulations was presented by Bournaud et 
al. (2004), who investigated tidal tails in interacting galax- 
ies, also with respect to the signatures they leave in the 
position- velocity diagram of the galaxy. They found that, 
depending on the viewing angle, the sign of the velocity 
gradient can change in the presence of tidal tails. 



This paper is organised as follows. In Sect. 2 and 3, we 
describe the simulations and the way we extract realistic 
rotation curves from the numerical data. Sect. 4 explains 
the interaction geometries used for the present work. In 
Sect. 5 the results for different interaction scenarios and 
their dependence on the viewing angle are presented, fol- 
lowed by a short discussion of some of the results' impli- 
cations for observational properties of cluster galaxies in 
Sect. 6. We conclude with a summary of the main conclu- 
sions in Sect. 7. 

2. Simulations 

In this work we use some of the interacting systems pre- 
sented in Kapferer et al. (2005). The simulations were 
carried out with the N-body/SPH code GADGET-2 de- 
veloped by V. Springel (see Springel 2005 for details). 
In this code the gas of the galaxies is treated hydrody- 
namically, including cooling, star formation, stellar feed- 
back, and galactic winds (Springel & Hernquist, 2003). 
The collisionless dynamics of the dark matter and the stel- 
lar component is calculated using an N-body technique. 
The combined N-body/SPH simulation then calculates 5 
Gyr of evolution. For every time step, we know the ve- 
locity of each particle and can extract realistic rotation 
curves. Other authors have used purely N-body techniques 
to investigate the kinematics of the modelled galaxies (e.g. 
Barton et al. il999t Heyl et al. 1996). Although the over- 
all kinematics of the disc is dominated by the darkmatter 
halo, disturbances of the gas have to be treated hydrody- 
namically, especially with respect to survival timescales of 
the distortions. The dynamics of the gas during an inter- 
action is significantly different from the dynamics of the 
collisionless stellar component. Thus, a purely collision- 
less approach is insufficient for describing the kinematics 
of the gas in detail. Note that also the star formation 
rate of the interacting system is usually higher than that 
of the isolated system (Kapferer et al., 2005). Therefore, 
stronger galactic winds can occur, which will, in general, 
decrease the rotational velocities (see paper I for a dis- 
cussion of this effect). Nevertheless, the results obtained 
by the above-cited groups are very interesting and in sev- 
eral points similar to ours, so we will use them frequently 
as reference in subsequent sections. lono et al. (2004) also 
used an N-body/SPH code to investigate position velocity 
diagrams of modelled galaxies, but focussed on signatures 
of radial gas flows. 

In Kapferer et al. (2006, hereafter paper I) we have 
shown that our models successfully reproduce rotation 
curves of isolated spiral galaxies. For completeness we re- 
peat the initial conditions of our model galaxies in Tabled 
These initial conditions were built according to Springel et 
al. (2005), which is based on the work of Mo et al. (1998). 
The two model galaxies were chosen such that they rep- 
resent a Milky Way type and a small spiral galaxy, where 
the mass ratio of the two galaxies is 8:1. 

Throughout the paper we adopt the standard ACDM 
cosmology with Oa = 0.7, fim = 0.3, and h=0.7. 
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Table 1. Properties of the initial conditions of the model 
galaxies. 



Properties 


Galaxy A 


Galaxy B 


circular velocity Vc^ 


160 


80 


disc mass fraction^ 


0.05 


0.05 


gas content in the disc^ 


0.25 


0.25 


disc thickness* 


0.02 


0.02 


total mass [Mq] 


1.33x10^^ h-^ 


1.67x10" h-^ 


disc scale length [kpc] 


4.51ft"^ 





circular velocity at 1200 in km/s 

fraction of disc particles (stars/gas) in units of halo mass 

relative content of gas in the disc 

thickness of the disc in units of radial scale length 




2.1. Advanced rotation curve extraction 

To follow the procedure applied in observations of distant, 
small, and faint galaxies as closely as possible we define as 
a first step a slit as shown in Fig. ^ Again, as in paper I, 
we choose large relative slit widths as they are used in ob- 
servations of intermediate and high redshift galaxies (e.g. 
Bohm et al. 120041 Nakamura et al. 12006(1 . As a next step 
we project all stellar and gas particles of the N-body/SPH 
simulations onto a cartesian, uniform grid, with a spacing 
corresponding to typical spatial resolutions of the observa- 
tions of intermediate redshift galaxies, for example, in the 
case of FORSl the pixel scale of 0.2", which corresponds 
to ^ 1.1 kpc for a redshift of z=0.4. Subsequently we 
derive a stellar blue band luminosity distribution J(x, z) 
(the coordinate system is chosen as indicated in Fig. ^ 
by assuming a stellar mass-to- light ratio of 1.2. This is 
the mean obtained using different star formation histories 
in the redshift range z = 0.5 — 1.4 in Dickinson et al. 
(2003). Then we calculate a luminosity- weighted velocity 
field Kv as: 

Vy,{x,z) = V{x,z) X I{x,z), (1) 

where V{x, z) is the projected, 2D velocity field of the 
gas. This weighting assures that brighter regions con- 
tribute more to the velocity shift. To simulate seeing ef- 
fects on our velocity-field measurements, a convolution 
with a Gaussian point spread function (PSF) was adopted 
on the luminosity-weighted velocity field with an FWHM 
according to typical observational conditions: 

Kv(x,^)^Kvc(a;,z). (2) 

The FWHM of the seeing in observations is typically of the 
order of 0.5"-1.2" (see e.g. Jager et al. l2()04|l . For our sim- 
ulations we adopt a mean value of 0.8" that corresponds, 
for example, to 4.3 kpc for a source at redshift z=0.4. 

Finally, we compute the mean velocity at a certain 
galactocentric radius r, by averaging over all grid points 
within the slit perpendicular to the spatial axis. Hence, 
the velocity extracted for the RC at radius r reads: 



Fig. 1. Model galaxy A (gas (blue) and stars (yellow) are 
shown) and a virtual slit for extracting a rotation curve. 
The slit width d is indicated by the red lines. A coordinate 
system as used in the text (see Sect. 2.1) is shown. 



Here a; = r is the apparent galactocentric radius along the 
spatial axis, and the sum runs over all grid points that 
lie within the defined slit. For all RCs presented in this 
work, we place the virtual slit along the major axis of the 
investigated galaxy. The normalisation with the convolved 
luminosity distribution Ic{x^ z) is necessary to ensure that 
the rotational velocity V{r) of a constant velocity field 
would not be changed by the weighting and the seeing. 
This approach allows a direct comparison with 2D obser- 
vations. In Fig. 13 we show the luminosity field derived 
from the stellar mass distribution, the unsmoothed and 
unweighted velocity field of the gas, and finally the ex- 
tracted RC for one of our modelled galaxies in an early 
stage of an interaction (unequal mass merger). 

We emphasise that the B-band mass-to-light ratio 
(M/L) is in general not constant for typical spiral galax- 
ies but a complex function of radius (e.g. Verheijen et al. 
2005). Therefore we extracted rotation curves by assuming 
different radial gradients in the mass-to-light ratio finding 
that the overall shape of the RC, which we investigate in 
this paper, is not sensitive to the specific choice of M/L. 
In Fig. Owe show RCs extracted from one of our modelled 
systems using two different assumptions for the mass-to- 
light ratio. In one case we used a constant mass-to-light 
ratio while in the other case we extracted the RC using a 
linearly decreasing M/L with a slope of 0.1 kpc~^. Some 
quantitative differences are visible, but the shape of the 
RC is not affected. 

2.2. The spatial alignment and impact parameters 

We follow the notation introduced by Due et al. (2000) 
to describe the interaction geometry. Parameter b corre- 
sponds to the minimum separation of the galaxies' trajec- 
tories, as if they were point masses on Keplerian orbits 
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Fig. 2. Top left: luminosity field derived from the stellar- 
mass distribution; top right: the unsmoothed and un- 
weighted velocity field of the gas. A virtual slit is indicated 
and a small part of the gaseous disc of the second galaxy is 
visible on the right hand side. Bottom left: the extracted 
RC for the modelled galaxy (dashed line) with la errors 
(solid lines). For this image series we use galaxy A in an 
early stage of an interaction (unequal mass merger). 



(see Fig. 0)). The two angles, & and <&, that define the 
spatial orientation of the disc are described in Fig. 01 

For the complete sample, we selected the alignments 
in such a way as to cover as many geometries as possible, 
including minor and major mergers and fly-bys (achieved 
by increasing the minimum separation). The appearance 
of intrinsic RC distortions will of course not only depend 



80 




-120 I ' • • 1 

-10 -5 5 10 

position [kpc] 

Fig. 3. Rotation curve extracted from one of our modelled 
systems along the major axis using two different assump- 
tions for the mass-to-light ratio. The dashed line corre- 
sponds to a constant mass-to-light ratio, while the solid 
line was extracted using a linear decreasing M/L with a 
slope of 0.1 kpc~^ 



on the interaction geometry but also on the angle under 
which the galaxy is observed (see Sect. 3.4). The interac- 
tion geometries we use for this work are listed in Table 13 
which represents a small sub-sample of the large sample 
introduced by Kapferer et al. (2005). 

3. Results 

In the following subsections we study the kinematic distor- 
tions caused by unequal mass mergers, equal mass merg- 
ers, and fly-bys and investigate the dependence of their 
appearance on the viewing angle. 
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Table 2. Interaction parameters for the simulations. C2: 
the two interacting model galaxies; C7: minimum separa- 
tion [kpc]; C8: initial relative velocities of the two galaxies 
[km/s]. 



Simulation 




$1 






H2 






1 


A R 

r\.-D 


u 


u 


n 
u 


n 
u 


c; 



zou 


2 


B-B 














5 


120 


3 


B-B 














50 


120 


4 


B-B 














25 


120 


5 


A-B 








180 





5 


250 


6 


A-B 














25 


250 


7 


A-B 














50 


250 


8 


A-B 








90 





5 


250 


9 


B-B 








90 





5 


120 


10 


B-B 








180 





5 


120 



3.1. Rotation curve distortions caused by unequal 
mass mergers 

To investigate the effects of unequal mass mergers on the 
RC of the massive partner in the interaction, we use sim- 
ulation 1 (see Table [SJ , which is a co-rotating coplanar 
collision of model galaxy A with model galaxy B. Figure 
|31 shows the Gaussian smoothed images of the luminosity 
field of the galaxies, derived from the stellar mass dis- 
tribution and with the assumption of a constant mass- 
to-light ratio (see Sect. 2.1). For the spatial resolution 
of the images, we chose 0.1"/pixel (0.54 kpc for z—Q.A), 
which corresponds to the optical resolution FWHM of the 
ACS/WFC instrument on the Hubble Space Telescope 
(HST). It is obvious from these images that minor mergers 
are not always easy to detect in direct imaging, especially 
for intermediate and high redshift galaxies. In the early 
stages of an interaction, it is often possible to resolve the 
galaxy pairs, so the interaction is clearly visible. Later 
on it is more difficult to distinguish between an isolated 
galaxy and a late stage of interaction. 

In Fig. El we show a sequence of images corresponding 
to different time steps in the merging process, where only 
the gas is shown. Note that the system is always viewed 
nearly face-on in this figure, while the RC is extracted 
from the system observed under an inclination angle of z = 
80° (edge-on is defined hy i — 90°). Chart (a) represents 
0.5 Gyr, (b) 0.6 Gyr, (c) 1 Gyr, and (d) 2 Gyr of evolution. 
In chart (a) in Fig. El the galaxy is shown shortly 100 
Myr) before the first encounter of the two gas discs, i.e. 
the moment, in which the two gas discs overlap completely 
for the first time. In the next chart (b) the two galaxies 
have their first encounter, whereas charts (c) and (d) show 
the galaxies 0.4 Gyr and 1.4 Gyr later, respectively. The 
massive partner (galaxy A) shows the typical features of 
merging events such as tidal tails and bridges. After 1.5 
Gyr the perturbation in the gas distribution is still visible 
in a pronounced spiral arm, as shown in Fig.O (d) on the 
right hand side. 



a 


b 


> 





c 


d 




• 







Fig. 5. Gaussian smoothed images of the luminosity field 
of the galaxies, derived from the stellar mass distribution 
and with the assumption of a constant mass-to-light ratio. 
The images (a) to (d) correspond to the images in Fig. |B1 
but observed under an inclination angle of 10°. The se- 
quence shows an unequal mass merger, where image (a) 
represents 0.5 Gyr, (b) 0.6 Gyr, (c) 1 Gyr, and (d) 2 Gyr 
of evolution. More details about the smoothing are listed 
in Sect. 2.1. 



In Fig. Ul the corresponding RCs are shown. The RCs 
were extracted as described in Sect. 2.1 for model galaxy 
A. Before the first encounter, the RC has an undisturbed 
shape. As galaxy B comes closer, galaxy A feels an accel- 
eration towards the common centre of mass. This acceler- 
ation affects the velocity field in such a way that relative 
velocities with respect to the centre of galaxy A are en- 
hanced towards the approaching galaxy B and weakened 
on the opposite side. In other words, the kinematic centre 
of the galaxy is shifted with respect to the geometrical cen- 
tre of the gaseous disc. The RC of the merging-phase (b) 
shows strong signatures of the first encounter in its asym- 
metric shape. The incorporated galaxy B with its velocity 
component along the line of sight causes the acceleration 
on the left hand side of the RC. The right hand side of 
the RC is not affected yet. This part of the RC is not 
distorted and does resemble the 'classical' rotation curve 
shape rising in the inner part and turning over to a fiat 
regime. The third RC corresponds to a later evolution- 
ary step (400 Myr after RC (a)), where the less massive 
galaxy has passed galaxy A and left behind a disturbed 
velocity field. Although the RC has become almost sym- 
metric again, it shows signs of interaction, especially the 
pronounced bumps in the outer parts (± ~ 10 kpc). The 
last RC (d) shows the system 1.4 Gyr after the first en- 
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counter; note that galaxy B is no longer visible in the 
direct vicinity. Indeed the system has almost relaxed and 
nearly no signs of the interaction are left. 

In order to quantify the degree of asymmetry of the RC 
shape, we follow an approach used in Dale et al. (201) 1 i. 
We measure the area between the kinematically folded ap- 
proaching and receding halves and divide it by the average 
area under the RC: 



Asymmetry = 



\\Vir)\-\V{-r)\\ 



\V{r)\ + \Vi-r)\ 



(4) 



Here cr(r) is the uncertainty of the rotational velocity 
at position r. We adopted a value of 20 km/s, which is typ- 
ically the minimum uncertainty in observations of distant 
galaxies due to the limited spectral resolution. In Table 01 
we list the asymmetry measure for different instances un- 
der different viewing angles for the unequal mass merger 
of simulation 1. The asymmetry of the RC increases with 
time as the two galaxies approach each other. After 0.5 
Gyr a significant asymmetry is present along all lines-of- 
sight. The asymmetry reaches a maximum after 0.7 Gyr 
of evolution and subsequently falls to the value it had be- 
fore the first encounter. The small drop in the asymmetry 
measure after 0.9 Gyr corresponds to the time between 
the first and second encounter. Note that the asymmetry 
strongly depends on the viewing angle at a certain instant. 
We will discuss this issue in Sect. 13.41 The small compan- 
ion galaxy B, on the contrary, already loses most of its gas 
during the first pass through galaxy A and finally merges 
completely with the larger galaxy; hence, the distortions 
in the velocity field are so severe that we cannot make a 
more detailed investigation. 

From both theoretical and observational studies, it is 
expected that kinematic disturbances fade within a few ro- 
tation cycles (<1 Gyr) and therefore reflect the recent in- 
teraction history (e.g. Rubin et al. 119991 Dale et al. 2001). 
Our results from this and the subsequent sections con- 
firm this statement numerically within the investigated 
interaction scenarios. The damping of the disturbances is 
listed quantitatively in Table 13 1.4 Gyr after the first 
encounter, i.e. after 2 Gyr of evolution, the asymmetry 
measure reaches a value close to the one after 0.1 Gyr of 
evolution, i.e. before the first encounter. 

In Fig. IHlthe RC for model galaxy A 1.4 Gyr after 
the first encounter, i.e. after 2 Gyr of evolution is com- 
pared to the RC of the isolated model galaxy A after 2 
Gyr of isolated evolution. The shape of the RC for model 
galaxy A resembles ~1 Gyr after the first encounter the 
one of an isolated system, so the properties derived from 
it would not differ significantly within given errors. The 
only significant difference is the slope of the inner part of 
the RC. For intermediate and high-redshift galaxies, small 
bumps as in the right section of the red line would not be 
observable. 



Table 3. Evolution of the asymmetry parameter^ for un- 
equal mass mergers (simulations 1, 5, and 8) under differ- 
ent lines-of-sight (angle $ is measured with respect to an 
arbitrarily chosen line-of-sight). 



Sim. 1: 




time 


100 


500 


700 


900 


1100 


2000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


9% 


54% 


103% 


20% 


100% 


12% 


45° 


10% 


44% 


120% 


58% 


103% 


17% 


90° 


7% 


21% 


112% 


67% 


89% 


8% 


135° 


8% 


52% 


45% 


46% 


51% 


7% 


Sim. 5: 




time 


100 


500 


700 


900 


1100 


2000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


10% 


50% 


48% 


20% 


45% 


6% 


45° 


10% 


11% 


51% 


13% 


12% 


18% 


90° 


12% 


70% 


38% 


15% 


33% 


12% 


135° 


9% 


61% 


12% 


9% 


31% 


17% 


Sim. 8: 




time 


100 


500 


700 


900 


1100 


2000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


11% 


21% 


122% 


44% 


19% 


13% 


45° 


9% 


93% 


83% 


38% 


26% 


17% 


90° 


13% 


74% 


40% 


13% 


28% 


20% 


135° 


10% 


110% 


38% 


10% 


37% 


7% 



m 



.asymmetry measure as defined by Dale et al. 2001, see Eq. 






Fig. 6. Sequence of images corresponding to different time 
steps in the merging process, with only the gas shown. The 
evolution of the system is calculated over 5 Gyr. Here (a) 
represents 0.5 Gyr, (b) 0.6 Gyr, (c) 1 Gyr, and (d) 2 Gyr 
of evolution. 



We study effects caused by collisions of counter- 
rotating discs using simulation 5 (see Table The shape 
of the RCs do differ from those of the co-rotating case, 
as different velocity fields are superposed. However, there 
are also velocity components due to the movement of 
galaxy B on its orbit that is approaching galaxy A. These 
are the same for both configurations and are dominating 
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-15 [kpc] 15 -15 [kpc] 15 



Fig. 7. Rotation curves for an unequal mass merger 
(galaxies A and B) for different time steps, correspond- 
ing to the time series shown in Fig. but observed under 
an inclination of i=80°. The red lines always indicate the 
centre of galaxy A to highlight asymmetries in the RCs. 



200 




_200i ' ' • • • • • 1 

-15 -10 -5 5 10 15 

[kpc] 

Fig. 8. Rotation curve of model galaxy A (blue, dashed 
line) 1.4 Gyr after the first encounter with model galaxy 
B, i.e. after 2 Gyr of evolution, as compared to the rotation 
curve of the isolated model galaxy A (red, solid line) after 
2 Gyr of isolated evolution. 



for the unequal mass merger, which we investigate here. 
Connected to the relative trajectories of the galaxies is also 
a transfer of angular momentum, leading to increasing ro- 
tational velocities (Heyl et al. 1996); therefore, the shape 
of the RC of galaxy A during the interaction shows only a 
small dependence on the rotational direction of galaxy B. 



3.2. Rotation curve distortions caused by equal mass 
mergers 

In a next step we investigate the effects of equal mass 
mergers on the RCs of the interacting system. For this 
purpose we use simulation 2 (see Table EJ, which is a co- 
rotating coplanar collision of two model galaxies B. Equal 
mass mergers massively disturb the velocity fields of the 
involved galaxies. Typically, the gaseous disc is almost 
completely destroyed in the merging process, and the or- 
bits of the stars are randomized, leading to the formation 
of an elliptical galaxy. This result has been shown in many 
previous studies (e.g. Barnes & Hernquist 1991; Cretton 
et al. 2001; Naab & Burkert 2003; Bournaud et al. 2005). 
We also find this formation of an elliptical galaxy in our 
simulations; however, as this issue was investigated com- 
prehensively by other authors, we do not discuss this in 
more detail in the present work. Instead, we chose two 
snapshots for which the gaseous disc is still observable 
and extracted the rotation curves for one of the interact- 
ing galaxies. The result is shown in Fig. |51 together with 
the corresponding interaction geometry. Note that in this 
figure the geometry is again illustrated as observed from 
nearly face-on, while the RC is extracted from the system 
observed under an inclination angle oi i = 80° , as in the 
previous section. Panel (a) shows the system after 1.1 Gyr 
and panel (b) after 2 Gyr of evolution. The first encounter 
of the two gas discs happens after 1.2 Gyr of evolution, 
i.e. 0.1 Gyr after the snapshot shown in panel (a). 

Before the first encounter, the RC of the investigated 
galaxy is only slightly disturbed but shows an asymmet- 
ric shape. As in the previous section, the kinematic centre 
of the galaxy is shifted with respect to the geometrical 
centre, i.e. the maximum of the surface mass density of 
the gaseous disc, due to the influence of the approaching 
galaxy. As the two discs pass through each other, tidal tails 
and the bridges are formed. These features can be seen in 
panel (b), where the distribution of the gas is shown. An 
RC extracted from the highlighted galaxy does not show a 
symmetric shape either. Again the merging process intro- 
duces a strong acceleration on the right hand side in the 
galaxy, leading to the formation of a tidal tail. Therefore, 
the relative velocities on the right hand side decrease with 
respect to the centre, while the left hand side shows higher 
relative velocities, leading to the asymmetric shape. The 
decrease in the rotational velocities in the outermost parts 
on the right hand side of the RC emerges as gas from the 
tidal tail is covered by the slit. The gas and the stars of 
the tail falling back onto the rotating disc have a small 
radial velocity component at this projected position, as il- 
lustrated in Fig.^l This possible change in the sign of the 
velocity gradients in the presence of tidal tails has already 
been found by Bournaud et al. (2004), who also points out 
the strong sensitivity of this result on the viewing angle. 
(We investigate this dependence on the viewing angle in 
Sect. 13.4] ) A significant decrease in the outer parts of the 
RC (see Fig. could indicate the presence of a massive 
tidal tail as consequence of a recent major merger event. 
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Fig. 9. Rotation curves for an equal mass merger for dif- 
ferent time steps corresponding to the interaction geome- 
try illustrated in the lower panel. The galaxy from which 
the RC was extracted is highlighted. Here panel (a) shows 
the system after 1.1 Gyr and panel (b) after 2 Gyr of 
evolution. 



However, this feature is not unique for tidal tails. As we 
will show in Sect. 13.41 the viewing angle of the observa- 
tion also plays an important role (see also Bournaud et al. 
2004). 

3.3. Rotation curve distortions caused by fly-bys 

In order to estimate the influence of fly-bys on the shape 
of RCs, we use simulation 3 (see Table which has a 
minimum separation of 50 kpc (at the first encounter). 
In Fig. ^2 the RC of the highlighted galaxy for different 
snapshots is presented in the upper panel. In the lower 
panel the geometry corresponding to the different time 
steps of the interaction is shown. Panel (a) represents 0.5 
Gyr, (b) 1.7 Gyr, (c) 2 Gyr, and (d) 3 Gyr of evolution. 
The change in the velocity field of the gaseous disc on 
short timescales is apparent in the right panel. The geom- 
etry (b) represents the nearest encounter, which is visible 
by a highly asymmetric shape. On the side pointing to- 
wards the companion, the RC does not flatten, due to the 
additional acceleration. On the opposite side, this feature 
is not present, instead the RC flattens. After 300 Myr 
(c), the velocity field does not show this asymmetry any- 
more, although the overall shape does not resemble an 
undisturbed RC. In the left panel we show the RC 0.5 



Fig. 10. Sketch to illustrate the decrease in the rotational 
velocities in the outer parts of the RCs that can be caused 
by gas from the tidal tail, which has small radial velocity 
components at this projected position and is covered by 
the slit. 



Gyr after the beginning of the simulation (a) and after 
2.5 Gyr of evolution (d). Note that the galaxies are spa- 
tially well-separated in both cases and that they would be 
therefore considered as isolated systems especially at high 
redshifts. The fly- by does not lead to persistent changes in 
the velocity field, and consequently the RCs do not show 
significant differences before and after the fly-by. The clos- 
est approach and the formation of observable distortions 
happens after ~ 2 Gyr of evolution (panels (b) and (c)). 
After 3 Gyr of evolution, i.e. 1 Gyr after the appearance of 
the flrst distortions, the RC no longer shows signs of inter- 
action. Thus, the velocity fleld again settles in this case to 
a relatively undisturbed equilibrium state after ^1 Gyr. 

Again we quantify the degree of asymmetry using Eq. 
Eland present the results in Tabled] We find an increas- 
ing value for the asymmetry with time as the two galaxies 
approach each other. After ^ 2 Gyr of evolution, the two 
galaxies come closest to each other. One Gyr later, the 
RC of the investigated galaxies shows a relatively undis- 
turbed form again under most viewing angles, which again 
quantifies the settling time of ^ 1 Gyr. In this interaction 
scenario we also find the strong dependence on the viewing 
angle in an interacting system. 

Decreasing the minimum separation to 25 kpc as in 
simulation 4 (see Table 0) leads to slightly stronger dis- 
tortions during the closest approach of the two galax- 
ies, both in the mass distribution and in the kinematics. 
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Fig. 11. Rotation curves for different snapshots of a fly- by 
of two model galaxies B. In the lower panel the alignment 
of the system and in the upper panel the corresponding 
rotation curves of the highlighted galaxy are shown. Here 
(a) represents 0.5 Gyr, (b) 1.7 Gyr, (c) 2 Gyr, and (d) 3 
Gyr of evolution. 



Nevertheless, the system settles back to an equilibrium 
state on approximately the same timescale as in the more 
distant fly- by discussed above (~ I Gyr). A fly- by of a less 
massive galaxy leads to lower asymmetries in the rotation 
curve (see simulation 7 in Table 0J. 

3.4. Dependence of the results on the viewing angle 

In an interacting system, two different velocity fields over- 
lap. Therefore, the shape of the RC, which is a ID rep- 
resentation of the full 3D velocity field, will significantly 
depend on the viewing angle. This fact is shown in Fig. 1121 
for the RCs extracted from the first encounter of model 
galaxies A and B seen nearly edge on from different view- 
ing angles. The geometry is illustrated by the image series 
on the left hand side of the figure. Again this image series 
shows the system face-on, while the RCs are extracted 
from the galaxy observed under an inclination angle of 
i = 80°. We varied the azimuthal angle in steps of 45°. 
The differences in the RCs are striking and are mainly 
caused by the velocity components of galaxy B on its tra- 
jectory and its projections. 

As we have shown in Sect. 13.21 the presence of a tidal 
tail can lead to a change in the sign of the velocity gradi- 
ent, which was also described by Bournaud et al. ^2.004 ). 
However, this feature depends on the viewing angle, as in 
long-slit spectra (and also in position-velocity diagrams 
obtained from Fabry-Perot intcrferomctry) we always look 
at only one (the linc-of-sight) component of a 3D vector 
field. If this component is smaller at a certain position, 
because, for example, the gas there has been decoupled 
from the general rotation of the disc by tidal interaction, 
a decrease in the RC can be found. In Fig. El we show 
the interacting system with pronounced tidal tails from 



different lines-of-sight. Note that even a flat RC is possi- 
ble in some configurations, although we found a change 
in the sign of the velocity gradient for most lines-of-sight. 
We also note that a prominent turn over in the RC is 
present in the third instance shown in Fig. El A part 
of the tidal tail is aligned with the line-of-sight. In their 
simulations for such configurations Bournaud et al. H2004|l 
always found a change in the sign of the velocity gradient 
in their position- velocity diagrams, consistent with our re- 
sults shown in Fig. El In addition to the line-of-sight, the 
time evolution of the tidal tail also plays an important role 
in the question of whether the turn over in the RC can be 
observed. Hibbard & Mihos 1)1995(1 showed that the kine- 
matics of the tidal material depends on the radius, as well 
as on the time evolution. A decreasing outer part in the 
RC at a certain instant and under a certain line-of-sight 
can be interpreted, as pointed out in Fig. El However, 
this feature is not unique, as we mentioned before and a 
tidal tail is just one possible explanation. 

As there is an infinite number of viewing angles to 
an infinite number of different interaction geometries, it 
is impossible to relate a specific observed RC distortion 
unambiguously to a certain interaction scenario without 
further information e.g. from high-resolution direct imag- 
ing. The specific shape of the RC will always depend on 
both the geometry of the interaction and the viewing an- 
gle. The situation is further complicated, as the dynami- 
cal evolution of the system is highly non-linear. However, 
these interaction-induced distortions have in common that 
the RCs become asymmetric, show significantly rising or 
falling profiles on the side pointing towards the other 
galaxy, and show pronounced bumps depending on the 
viewing angle (cf. Figs. 171 and I12|) . We found no observa- 
tional biases for isolated galaxies that could produce such 
pronounced features (see paper I for details). Most obser- 
vational effects change the RC of a galaxy in a symmet- 
ric way (e.g. relative large slit widths, inclination). Spiral 
arms typically cause fluctuations of a few tens of km/s 
that are superposed on the smooth rotation curve of the 
galaxy (see e.g. Sofue and Rubin, 2001 and paper I). Thus 
spiral arms can only account for smaller distortions but 
not for high asymmetries, as present in the present paper. 
Of course it is important to stress that galaxies in the uni- 
verse are never isolated. The dynamical history of a galaxy 
will always influence the kinematics of the different com- 
ponents of galaxies. Not all asymmetries are caused by 
clearly observable mergers/interactions (Emergreen 2005 
and references therein), although distinct asymmetries are 
often observed in interacting systems. Very recently dis- 
turbed isolated galaxies were investigated to search for in- 
dicators of interactions with dark galaxies, i.e. DM halos 
without any baryonic matter (Karachentsev et al. 2006). 

4. Implications of the results for the observational 
properties of cluster galaxies 

If the results, which were found in the previous sections 
for a restricted sample of interaction geometries, hold 
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Table 4. Evolution of the asymmetry parameter^ for fly-bys (simulations 3, 4, and 7) under different lines-of-sight 
(angle $ is measured with respect to an arbitrarily chosen line-of-sight). 



Sim. 3: 




time 


700 


900 


1100 


1300 


1500 


1700 


2000 


3000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


12% 


25% 


45% 


49% 


89% 


99% 


44% 


14% 


45° 


4% 


19% 


35% 


59% 


35% 


30% 


19% 


9% 


90° 


15% 


21% 


15% 


62% 


16% 


16% 


17% 


4% 


135° 


20% 


18% 


21% 


19% 


97% 


2% 


4% 


23% 


Sim. 4: 




time 


700 


900 


1100 


1300 


1500 


1700 


2000 


3000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


16% 


17% 


17% 


16% 


37% 


20% 


15% 


8% 


45° 


25% 


16% 


31% 


19% 


16% 


53% 


16% 


5% 


90° 


26% 


17% 


48% 


22% 


40% 


20% 


21% 


10% 


135° 


30% 


10% 


21% 


24% 


36% 


90% 


16% 


20% 


Sim. 7: 




time 


700 


900 


1100 


1300 


1500 


1700 


2000 


3000 


$ 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 


0° 


8% 


15% 


17% 


6% 


9% 


12% 


13% 


5% 


45° 


20% 


8% 


14% 


13% 


8% 


15% 


9% 


10% 


90° 


21% 


13% 


16% 


20% 


13% 


14% 


7% 


12% 


135° 


17% 


14% 


12% 


21% 


15% 


15% 


6% 


10% 



^....asymmetry measure as defined by Dale et al. 2001, see Eq.0] 
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Fig. 12. Rotation curves for the first encounter of model 
galaxy A and B seen nearly edge on (i=80°) from different 
viewing angles. The geometry is illustrated by the image 
series on the left hand side. The snapshot corresponds to 
panel (b) in Fig. and therefore shows the system after 
0.6 Gyr of evolution. 




line of 
sight 



Fig. 13. Rotation curves for model galaxy A after an equal 
mass merger from different viewing angles. The geometry 
is illustrated by the image series on the left hand side. The 
snapshot corresponds to panel (b) in Fig. Inland therefore 
shows the system after 2 Gyr of evolution. The change in 
the sign of the velocity gradient does not appear for all 
viewing angles. 



more generally, then they have an important implication 
for observational properties of cluster galaxy populations. 
First of all, stable regular discs show the classical RC, ris- 
ing in the inner part and turning over to a flat regime. 
Interactions as mergers or fly-bys introduce instabilities 
to the gas, which lead to an enhanced star formation rate 
(e.g. Kapferer et al., 2005), provided the galaxy has not 
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lost large fractions of its gas e.g. by ram-pressure strip- 
ping. As shown in the previous sections, these interac- 
tions also cause kinematic distortions. After the interac- 
tion, both the star formation and the kinematic distor- 
tions again settle back into a stable equilibrium state. 
Our results therefore suggest that enhanced star forma- 
tion in apparently isolated galaxies can be triggered by 
non-equal mass mergers. In direct imaging, especially for 
high-redshift galaxies, those interactions are hardly vis- 
ible. Additional investigation of the velocity field would 
help to reveal the dynamical state of the system. This im- 
plies that distortions in RCs can give evidence of the star 
formation of a galaxy and vice versa. As consequences of 
interactions being even longer observable in the velocity 
field than in direct imaging, these investigations can help 
to extract clues for the star formation history of a galaxy. 

We now discuss briefly to what extent observations of 
distant cluster galaxies can be interpreted by simulations 
as presented here. To that purpose, we show thumbnails 
and rotation curves of four galaxies in Fig. E| as exam- 
ples. They are members of the cluster MS 0451.5-0305 at 
a mean redshift of z = 0.53, which was recently imaged 
by Ellis et al. using the ACS camera onboard HST in the 
F814W filter. The VLT/FORS spectra were obtained as 
part of a larger program to investigate galaxy transforma- 
tions as described by Ziegler et al. ( 2003I) (see also Jager et 
al.. l2()()4;i . Tilted MXU shtlets were placed along the pho- 
tometric major axis as determined on groundhased images 
(VLT/FORS) since no full coverage of the cluster field by 
HST was available yet. The slitlets with 0.9" width are 
actually much longer than indicated on the thumbnails. 
The position- velocity diagrams (still unpublished) present 
the measured shifts of emission lines using the method in- 
troduced by Bohm et al. H2004(l and Jager et al. H2004|l 
so that the displayed velocities are not corrected for the 
galaxy inclination and the effects caused by the observa- 
tions. Note, however, that our Vmax measurements fully 
accounted for geometrical (e.g. inclination) and observa- 
tional effects (like seeing). 

Two of the galaxies have rather regular structures and 
rotation curves, while the other two spirals are clearly dis- 
torted. They exhibit features that are very similar to what 
is seen in our simulations. 

5. Summary and conclusions 

In this work we have investigated the rotation curves 
of interacting spiral galaxies using N-body/SPH simula- 
tions with included prescriptions for cooling, star forma- 
tion, and stellar feedback. We focussed on the questions 
of which distortions of the RCs can arise from mergers (in 
unequal mass mergers the mass ratio of our galaxies is 8:1) 
and fly-bys (25 - 50 kpc), how long these disturbances sur- 
vive, and what is their dependence on the viewing angles. 
To summarise and conclude: 

— Galaxy-galaxy mergers significantly disturb the veloc- 
ity fields and hence the RCs of the interacting galaxies. 



These distortions commonly introduce asymmetries in 
the RCs with significantly rising or falling profiles on 
the side pointing towards the companion galaxy and 
pronounced bumps. These distortions are also clearly 
observable for distant galaxies. 

— The specific shape of the RC depends on the mass ratio 
of the interacting galaxies, the geometry of the inter- 
action, and the viewing angle. Therefore, without any 
additional information on the system, it is impossible 
to relate a specific, observed RC distortion unambigu- 
ously to a certain interaction scenario. High-resolution 
direct imaging and 2D velocity fields allow further in- 
sight into the dynamical state of the system. 

— We used a quantitative measure for the RC asymme- 
try and find for given interaction scenarios that the 
disturbed velocity field settles relatively fast (~1 Gyr) 
to a stable undistorted equilibrium state after an un- 
equal mass merger or a fiy-by. This numerically con- 
firms results from both previous theoretical and ob- 
servational studies. Therefore, the shapes of the RCs 
before and after such an interaction are very similar, 
especially as there are no persistent distortions visible 
in the RCs. Investigations of the number of interac- 
tions as a function of redshift including velocity fields 
could be strongly biased by this effect. 

— Tidal tails, commonly formed by equal mass mergers, 
can be visible in the rotation curves as a decrease in 
the rotational velocities, provided the gas with a small 
radial velocity component lies within the slit. However, 
we showed that this change in the sign of the velocity 
gradient is not unique, and a tidal tail is just one of the 
possible explanations. Additionally, the viewing angle 
on the system plays an important role in the question 
of whether the turn over in the RC can be observed. 

— Two effects contribute to the distortions in the veloc- 
ity field of an interacting system, starting with the 
superposition of the two individual velocity fields that 
are modified by the disturbed gravitational potential. 
Also, the velocity components from the trajectories of 
the galaxies add to the resulting velocity field, and the 
relative orbits cause an angular momentum transfer 
(see also Heyl et al. 1996). 
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Fig. 14. ACS images taken in the F814 filter and rotation curves from the VLT/FORS spectroscopy of four galaxies 
in the distant cluster MS 0451.5-0305 at z=0.53. Tentatively, we find evidence for, e.g., tidal streams (upper panel 
left, cf. Fig. 8) or a recent close fly-by/merger (lower panel right, cf. Fig. 10). Note that the rotation curves give the 
observed values of the rotational velocity Vrot; i.e. they are not corrected for inclination or offset angles (S) between 
the slit directions and the galaxies' major axes. The indicated rectangles correspond to the extensions of the derived 
RCs. 
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